Clusters and Fluctuations at Mean-Field Critical Points and Spinodals 
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We show that the structure of the fluctuations close to spinodals and mean-field critical points is 
qualitatively different than the structure close to non-mean-field critical points. This difference has 
important implications for many areas including the formation of glasses in supercooled liquids. In 
particular, the divergence of the measured static structure function in near-mean-field systems close 
to the glass transition is suppressed relative to the mean-field prediction in systems for which a 
spatial symmetry is broken. 
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The structure of the fluctuations near critical points 
and spinodals and their relation to the behavior of quan- 
tities observed in experiments and simulations is impor- 
tant for understanding the properties of many materials. 
For critical phenomena in non-mean-field Ising systems 
this relation was found by mapping the thermal critical 
point onto a percolation transition [QJ|] . In these systems 
the properties of the clusters at the percolation threshold 
are identical to those of the fluctuations at the thermal 
critical point. In particular, the mean cluster diameter is 
the correlation length, the density of the spanning cluster 
scales as the order parameter, and the mean number of 
sites in the spanning cluster is the susceptibility [^],[| . 

For mean-field Ising models there is a line of spinodal 
critical points as well as the usual critical point. These 
mean-field thermal singularities also can be mapped onto 
percolation transitions Q , but the relation between per- 
colation clusters and critical fluctuations is qualitatively 
different. This difference has important consequences for 
supercooled mean-field liquids, which also exhibit spin- 
odals ||, and for "near- mean- field" systems, which are 
characterized by long, but finite range interactions and 
which exhibit many of the characteristics of mean-field 
systems including pseudospinodals. Many dense systems 
with short-range interactions exhibit near-mean-field be- 
havior under certain circumstances [§]|],@,|9|] ■ 

In the following we give scaling arguments that relate 
the structure of clusters at mean-field critical points and 
spinodals to measurements of thermal quantities such as 
the static structure function S(k) and discuss the impli- 
cations of these results for understanding mean-field and 
near-mean-field supercooled liquids. One motivation for 
this analysis is that the behavior of a supercooled two- 
component, two-dimensional (d — 2) Lennard- Jones sys- 
tem has been interpreted as due to the influence of a 
pseudospinodal (relative to the stable solid) |?],^| . A con- 
sequence of this pseudospinodal is that the first peak of 
S(k) is predicted to exhibit a power-law divergence, a 
prediction consistent with simulations in d = 2 . How- 
ever, experiments [nOl and simulations pdl] have failed to 



find similar behavior for S(k) in d = 3. This failure is 
puzzling because d — 3 Lennard-Jones systems should 
be more mean-field-like than d — 2 systems. This work 
addresses this apparent contradiction. 

We first review how the cluster structure relates to 
thermal critical phenomena in non-mean-field systems. 
In mapping the percolation transition onto the thermal 
critical point, the quantity that is isomorphic to the sin- 
gular part of the free energy is the mean number of clus- 
ters in a correlation region of volume £ d , where £ is the 
correlation length H . The singular part of the free energy 



scales as £ d e 2 



- dv , where e = (T - T c )/T c , and 



a and v are the specific heat and correlation length expo- 
nents respectively [0|. Because hyperscaling is valid, we 
have dv = 2 — a, and hence the mean number of clusters 
in a volume £ d is order unity B . 

To obtain the scaling of the isothermal compressibility 
or susceptibility xt, we use the relation of \t to the 
spatial integral of r(r), the order parameter correlation 
function [^3|. We take T(r) to be the square of the order 
parameter over a volume £ d and zero outside, and obtain 
Xt ~ e 2/5 £ d = e 2 - a - d "-~< using a + 2/3 | 7 = 2@. 
Hyperscaling gives 2 — a — dv = and hence xt ~ e 7 as 
expected. This simple argument is well known, but fails 
in mean-field systems where hyperscaling is not valid. 

In contrast to systems that obey hyperscaling, we need 
to distinguish between critical phenomena fluctuations 
and clusters in the mean- field limit. To discuss the latter 



we consider a weak long-range (Kac) interaction JL4 15 
of the form ^u^r), where u{pfr) is integrable, a function 
only of the distance between spins, and ferromagnetic. If 
the system size is first taken to infinity and then the range 
of interaction R = 7 _1 — > oo, we obtain the Curie- Weiss 
description of the thermodynamics and the Ornstein- 
Zernicke form for T(r) (ljjl^]. The exponents for the 
mean-field critical point are v = 1/2, a = 0,/3 = 1/2 and 
7=1 |Q . The singular part of the free energy at the 
mean-field critical point scales as 



A/ 



R 



d £ 2 



-d/2 



(1) 



1 



The quantity R d appears because all lengths are in units 
of the interaction range. 

From the percolation mapping Eq. ([!]) implies 

that the mean number of clusters in a volume £f is 
jj>d £ 2-d/2 yy e can estimate the magnitude of this num- 
ber from the Ginsburg criterion |l6) , which states that 
a system is well approximated by mean-field theory if 
XT/</> 2 £ d << 1, where ~ e 13 is the order parameter. If 
we use mean- field values for the critical exponents, we 
obtain the condition R d e 2 - d / 2 >> 1 f or mean-field the- 
ory to be valid. This inequality, together with Eq. (|l|) 
and the percolation mapping, implies that the number 
of clusters in a volume £ d is much greater than unity for 
mean-field systems. What does this result imply for the 
nature of clusters in mean-field systems? 

In non-mean-field systems which have one cluster in a 
volume £ d , the free energy cost of this cluster is ~ £ d g 2 ~". 
As discussed above, this cost is order unity in non-mean- 
field systems. In contrast, because A/ >> 1 in mean- 
field systems (see Eq. (Q)), e 2 ~" is not the free energy 

density of one cluster in mean-field. Rather, there are 
R d e 2-d/2 clusters 

in a volume £ d , and each cluster has a 
free energy cost of order unity. This free energy implies 
that the probability of a cluster is order unity. 

It also is clear that e^ 3 is not the density of one cluster, 
because if it were, the spin density in a volume of order 
£ d would be « R d e 2 ~ d ' 2 e^. For a fixed value of e, R can 
be made arbitrarily large and hence this density can be 
arbitrarily large, an absurd result. This argument and 
the one for A/ implies that e" is the density of all the 
spins in a volume £ d regardless of the cluster to which 
they belong. It also implies that the density of spins in 
one cluster near the mean-field critical point Wq is 



F V2 



R d £ 2-d/2 ' 



(2) 



This prediction has been verified numerically |18 
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Are these clusters the critical phenomena fluctuations? 
To answer this question we calculate the density of mean- 
field fluctuations. The partition function for the Gaus- 
sian approximation of the </> 4 model is given by 13 



Z=jS(f> exp{-/3 J dxR d [(\7<j)(x)) 2 + £0 2 (f)]}, 



(3) 



where (3 = l/ksT and ks is Boltzmann's constant. Be- 
cause we are interested only in the scaling properties of 
the fluctuations, we take the order parameter 4>(x) to be 
a constant 4> over a volume £ d and zero outside, and ob- 
tain Z ~ J5(j) e~P R e 1 $ . If we integrate <f> until the 
argument of the exponent becomes order unity, we find 
that the density of the critical phenomena fluctuations 
scales as 



,1/2 



( R d e 2-d/2y/2 



(4) 



That is, the structures with the density in Eq. (Q) have 
a free energy cost of one. As in the non-mean-field case, 
we expect these objects to be the critical fluctuations. 

From Eqs. (§) and (g) we see that the critical phenom- 
ena fluctuations are not the clusters, but are considerably 
denser in the mean-field limit R d e 2 - d / 2 >> 1. The struc- 
ture of the critical phenomena fluctuations in mean-field 
systems is that the "vacuum" is not featureless but con- 
tains a very large number (R d e 2 ~ d / 2 in a volume £ d ) of 
clusters. At the critical point (magnetic field ft, = in 
Ising models), the mean number of up and down clus- 
ters is equal giving a zero mean magnetization. Because 
these clusters are independent Jf|,^|, the fluctuations in 
the number of clusters is order (R d e 2 ~ d ^ 2 ) 1 ^ 2 . These fluc- 
tuations make up the zero magnetization background and 
are the critical phenomena fluctuations. If we use Eq. (Q), 
Xt - - eR d £- d / 2 /R d e 2 - d / 2 = e" 1 as expected. 

These arguments can be extended to spinodals in 
mean-field Ising models. If the spinodal is approached by 
varying h and keeping T fixed, the density of the clusters 
is Ah 1 / 2 /R d Ah 3 / 2 - d / i with Ah=h s -h^. The nu- 
merator represents the magnetic field scaling of the order 
parameter. The denominator is very large in mean-field 
and is the Ginsburg criterion for spinodals . As at the 
critical point, the mean number of up and down clusters 
is equal when the infinite cluster, which is related to the 
metastable magnetization, is subtracted. All scaling ar- 
guments for mean-field critical points apply at spinodals 
with appropriate changes in the values of the exponents 
(u= 1/4, a = 1/2,(3= 1/2, and 7 = 1/2 §§). 

We can determine the decay time of the clusters and 
the critical phenomena fluctuations by constructing an 
action from the linearized Langevin equation describing 
the dynamics of an Ising model with long-range inter- 
actions [^2). If we assume a random Gaussian noise, the 
probability measure for the order parameter cf>(x, t) is |23fl 



exp< — /3 / dxdt 



,dmt) 2 + {MR d e 2-d /2) 2 

at 



{-V 2 4>(x,t) + e^{x,t)} 2 + H@,rl>) 



(5) 



where M is a (constant) mobility, 



= e 1 / 2 ^, and 

#(V>,V0 = + MR d e 2 - d / 2 (-V 2 + e})^{x,t). 

The variables ^(x, t) and ip(%j t) obey a Grassmann alge- 
bra J22|;^3| and can be used to convert the average over 
the Gaussian noise to an average over the order param- 
eter p2|] . In the mean- field limit, the Langevin equa- 
tion describing the dynamics is linear, and the Grass- 
mann fields and the order parameter are independent 
p3[ . Hence, the Grassmann fields are irrelevant to a cal- 
culation of averages of functions of <fi. Using the action 
in Eq. (^|) , we can average a function of <fi by functionally 
integrating <fi up to its value at the critical point where we 
expect that the argument of the exponential in Eq. (|^) 
is order one. Because all of the terms in the action in 
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Eq. (||) are real and positive, each term in the exponen- 
tial must be order one, which implies that et; d <fi 2 /t ~ 1. 
Hence, the time scale that an object lives depends on cf). 

For critical fluctuations, ~ {R d t 2 ~ d / 2 )~ 1 / 2 , which 
leads to t ~ e -1 , the scaling for the decorrelation time 
in a mean-field system with a non-conserved order pa- 
rameter. For clusters, <f> c \ ~ (R d e 2 ~ d / 2 )~ 1 , which leads 
to Td ~ e^ 1 /R d e 2 ~ d ' 2 . We stress that these arguments 
are valid only in the mean-field limit and are a good ap- 
proximation in the near-mean-field case. 

These considerations imply that the clusters have a 
lifetime that is considerably shorter than the lifetime of 
the critical phenomena fluctuations. In the mean-field 
limit, R d e 2 ~ d l 2 — + oo, the lifetime of the clusters is zero. 

We now apply these ideas to a mean-field model of a 
supercooled fluid. We consider the step potential u(7 r ) = 
for 7r > 1 and u(~/r) = 1 for jr < 1 J24L For 7^0, 
S(k) = 1/[1 + p/3u(k)] in the fluid phase [5[, where p is 
the density and u(k) is the Fourier transform of -f d u(x) 
with x and k~ l scaled by 7. This mean- field form of 
S(k) indicates that the system is unstable for p/3 such 
that 1 + pj3 s u(ko) — 0, where u(ko) < is the minimum 
of u{k). As T -> T s = (k B f3 s )-\ S{k ) ~ e" 1 @, where 
e = \T — T s \/T s . The divergence of S(ko) is unchanged if 
a short-ranged reference potential is added. 

The scaling laws for the spinodal can be rewritten with 
e as the scaling variable |2f| with v = 1/2, 2 — a = 3, 
[3 = 1, and 7 = 1 rather than the values quoted previ- 
ously for the Ah scaling field EG] . The number of clusters 
in a volume £ d scales as R d e 3 -^/ 2 ( the density of clusters 
as e/R d e 3 ~ d / 2 , and the density of critical phenomena fluc- 
tuations as e/(i? d e 3_d / 2 ) 1 / 2 . The Ginsburg criterion is 
A = Ji d e 3-d/2 >> 1 anc j time scale for the critical 



and e~ 1 /R d e 



d 3-d/2 



fluctuations and the clusters is e~ 
respectively. 

The primary difference between the Ising and fluid 
spinodals is that S(k) diverges at fco 7^ for the fluid 
0. The Fourier transform of S(k) yields a correlation 
function T(r) that scales as r~ 1 e~ r ^ sin k$r in d = 3, 
implying that the critical fluctuations near the spinodal 
have a characteristic length k^ 1 « £ on which there is a 
periodic spatial variation. Since the critical fluctuations 
are an incoherent superposition of A overlapping clusters, 
the spatial symmetry breaking reflected in T(r) and S(k) 
also occurs for the clusters. We can show from an analysis 
of the Langevin equation that the clusters have a trian- 
gular structure in d = 2 and a bec or layered triangular 
structure in d = 3 ^5|. (Nuclcation near the spinodal is 
a coalescence of clusters |2l[] and the nucleation droplets 
have the symmetries discussed above p7|,p[.) 

To understand the consequences of our interpretation 
of mean- field and ncar-mean-field fluctuations, we obtain 
the T-dependence of S(ko) for a fluid using scaling argu- 
ments similar to the ones used above. The structure func- 
tion is related to an integral over the particles in a clus- 
ter times the probability that two particles belong to the 



same cluster within a critical fluctuation. The integral in- 
volves a phase factor ^ . exp[iko,j ■ Tj], where \ko,j\ = 
and the index j labels the basis reciprocal lattice direc- 
tions for the indicated symmetry. This phase factor is 
multiplied by the probability that a site belongs to a 
critical phenomena fluctuation, e/(i? d e 3_d / 2 ) 1 / 2 , times 
the probability that the second particle belongs to the 
same cluster, e/(R d e 3 ~ d / 2 ), times the number of clusters, 
(i? d e 3-d / 2 ) 1 / 2 . In addition, we average over the angles 
corresponding to the random orientations of the clusters. 
These considerations yield the scaling form: 



S(k ) 



dk df f R d e 3-d/2 ( 



ftd e 3-d/2 



£~i? 



l-d -(3-d)/2 



(6) 



The spatial integral is over a region the size of £ d and 
dk' denotes an integral over angles. 

Eq. (|)) predicts that S(ko) ~ with 7 = 1 in d = 1, 
7 = 1/2 in d = 2, and 7 = in d > 3 in contrast to 
7=1 for all d as predicted by mean-field theory ||. We 
stress that the weakening of the divergence of S(ko) in 
d = 2 and its supression for d > 2 is a consequence of 
the quenched periodic structure of the clusters. If this 
structure is modified by finite size effects or defects, then 
the suppression might not be as strong. 

The difference between the two calculations for the 
scaling behavior of S(ko) is the limiting procedure. If 
the limit A — > 00 ia taken before the calculation of S(ko), 
the lifetime of the clusters is zero in comparison to a mea- 
surement time. In this limit the clusters appear rotation- 
ally symmetric and a calculation similar to Eq. (^) would 
yield S(k ) ~ for all d, the same result as in Ref. ||. 
However, if we assume A to be arbitrarily large but fi- 
nite, the lifetime of the clusters is nonzero and the mea- 
surement time is smaller than the cluster lifetime. This 
assumption is consistent with the way measurements are 
done in experiments and simulations. 

It is difficult to estimate critical exponents by fitting 
data directly to the suggested asymptotic form. However, 
because the spinodal is well defined only in the mean- 
field limit and simulations can be done only for finite R, 
we must estimate 7 by approaching the pseudospinodal. 
But we cannot approach it too closely because we will 
reach the Becker-Doring limit where the system nucleates 
very quickly. Because the Metropolis algorithm becomes 
very inefficient in near-mean-field because almost all sin- 
gle particle moves are accepted, we included a small hard 
core of diameter a to decrease the acceptance probability. 

Our Monte Carlo results (see Fig. 1) in d — 1 were 
fit to AT^e-^ with 7 = 1.0 ± 0.05 and T s « 0.83. The 
T 7 term accounts in part for corrections to scaling which 
typically depend on d. The errors were estimated from 
different possible fits. The d — 2 results were fit to Ae~^ 
with 7 = 0.4 ± 0.05 and T s « 0.81. The quality of the 
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FIG. 1. The T-dependence of S{ko) in d = 1, 2, and 3 for 
the step potential discussed in the text. The d = 1 results 
are for N = 10 4 particles, interaction range R = 48, hard 
core diameter a = 0.05, and density p = 1.95. The fit shown 
for S(&o) is to Aft^ with 7 = 1.02 and T s = 0.82. For 
comparison, T 3 = 0.85 in the R — > 00 limit. The d = 2 results 
are for N = 10 3 , i? = 6, cr = 0.15, and p = 1.95. The fit for 
S(ko) is Ae~^ with 7 = 0.44 and T s — 0.81. For comparison, 
T 3 — 0.80 for R —* 00. The number of Monte Carlo steps per 
particle (MCS) for each value of T is 10 5 in d = 1 and 2. The 
d = 3 results are for AT = 10 3 , R = 3, cr = 0.45, and p = 1.95 
and fit equally well to Ae -7 with 7 = 0.16 and T s — 0.71 and 
to A log e + B with T s = 0.86; T s = 0.70 for i? -» 00. The two 
fits are indistinguishable in the figure. The number of MCS 
for each data point is between 2 x 10 5 and 4 x 10 s . 



fits is not as good as in d = 1, but the estimate for 7 
is consistent with our prediction of 7 = 0.5. The d = 3 
results are consistent with Ae^ with 7 = 0.16±0.02 and 
T s w 0.71 and with A log e + 5 with T s w 0.86. Details of 
the simulations will be discussed in a longer paper p5| . 

In summary, our theoretical and simulation results im- 
ply that the structure of the fluctuations in mean-field 
and near-mean-field systems differs qualitatively from 
that of non-mean-field systems. This structure leads to 
a suppression of the divergence of the measured static 
structure function near a pseudospinodal relative to the 
mean-field prediction in systems for which a spatial sym- 
metry is broken. The dependence on d of this suppression 
is such that there is no apparent divergence for d > 3 sub- 
ject to logarthmic corrections. In these systems there is 
a growing correlation length as the pseudospinodal is ap- 
proached that cannot be obtained from a measurement 
of S(k). We note that a divergent dynamical length has 
been found above the glass transition in the spherical p 
spin model p9| ]. This work, coupled with our results, 
raises the question of how the spinodal influences the dy- 
namical length. Our predictions have important impli- 
cations for the understanding of processes such as nucle- 
ation and glass formation in supercooled fluids. In par- 



ticular, we expect that the existence of clusters should 
help us understand the universal scaling behavior found 
for the dielectric response of organic glass formers . 
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